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ABSTRACT 

We investigate the nature of the newly discovered Ultra Faint dwarf spheroidal galax- 
ies (UF dSphs) in a general cosmological context simultaneously accounting for various 
"classical" dSphs and Milky Way properties including their Metallicity Distribution 
Function (MDF). To this aim we extend the merger tree approach previously devel- 
oped to include the presence of star-forming minihaloes, and a heuristic prescription 
for radiative feedback. The model successfully reproduces both the observed [Fe/H]- 
Luminosity relation and the mean MDF of UFs. In this picture UFs are the oldest, 
most dark matter-dominated [M/L > 100) dSphs with a total mass M = 10 ''^^Mq; 
they are leftovers of H2-cooling minihaloes formed at z > 8.5, i.e. before reionization. 
Their MDF is broader (because of a more prolonged SF) and shifted towards lower 
[Fe/H] (as a result of a lower gas metallicity at the time of formation) than that of 
classical dSphs. These systems are very ineffectively star-forming, turning into stars 
by z = only < 3% of the potentially available baryons. We provide a useful fit for 
the star formation efficiency of dSphs. 
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1 BACKGROUND 

. The dwarf spheroidal galaxy (dSph) population observed in 
' the Milky Way (MW) system provides crucial insights into 
\ cosmic structure formation. Nevertheless, the origin and evo- 
• lution of dSphs has so far escaped any comprehensive the- 
■ oretical interpretation and the picture is now further com- 
, plicated by the discovery of a new class of dwarf satellite 
galaxies: the Ultra Faint dSphs (UFs). 

UFs are the least luminous galaxies known, with a to- 
tal luminosity L « 1O"^~^L0; spectroscopic follow-up has re- 
vealed that they are highly dark matter dominated systems 
M/L > 100 (Simon & Geha 2007, Geha et al. 2008). Their 
average iron-abundance is ([Fe/H]) < —2 (Kirby et al. 2008) 
i. e. they represent the most metal-poor stellar systems ever 
known; although more data are required to solidly constrain 
stellar populations, at the moment all of them appear to be 
dominated by an old stellar population (Walsh, Willman & 
Jerjen 2008), with the only exception of LeoT (de Jong et al. 
2008). In addition, UFs are relatively common in the MW 
system, representing more than 50% of the total number of 
dSph companions discovered so far (A^ ~ 20); hence, they 
are not peculiar objects. When and how did UFs form? 

The [Fe/H] -Luminosity relation derived for UFs (Kirby 
et al. 2008) constitutes an extension towards lower metallic- 
ity of that of "classical" dSphs (L > lO^L©). Such a contin- 
uous trend seems to exclude that processes (e.g. tidal strip- 



ping) different from those shaping the relation for classical 
dSphs become dominant in these objects. However, while 
dSphs and UFs together span more than four orders of mag- 
nitude in luminosity, their total mass is roughly the same 
M ~ IO^Mq within the innermost 300 pc (Strigari et al. 
2008, Li et al. 2008). Why is then star formation so ineffi- 
cient in UF satellites? 

The puzzle is made even more intriguing by the obser- 
vation of metal-poor stars in UFs (Kirby et al. 2008) which 
have revealed the existence of a [Fe/H]< —3 stellar popula- 
tion, in contrast with the classical dSphs which are lacking 
extremely iron-poor stars (Helmi et al. 2006). Does this re- 
flect a different origin of UF and classical dSphs? In this Let- 
ter we extend previous works (Salvadori, Schneider, Ferrara 
2007 [SSF07]; Salvadori, Ferrara, Schneider 2008 [SFS08]) 
to investigate the origin of UF galaxies in a cosmological 
context. 



2 MODEL FEATURES 

In this Section we will briefly describe the basic features 
of the model implemented in the code GAMETE (GAlaxy 
MErger Tree & Evolution) which allows us to reconstruct 
the stellar population history and the chemical enrichment 
of the MW along its hierarchical merger tree, following at the 
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same time the formation and evolution of dSph satellites*. 
In such a model dSphs form out of the MW environment, or 
Galactic Medium (GM), whose metal enrichment strongly 
depends on the star formation (SF) history and mechanical 
feedback along the hierarchical tree. The model successfully 
reproduces the global properties of the MW and the Metal- 
licity Distribution Function (MDF) of Galactic halo stars, 
along with the observational properties of prototype dSphs, 
as for example Sculptor (SFS08). 

The merger history of the MW is reconstructed start- 
ing from z = 20 via a Monte Carlo algorithm based on the 
extended Press-Schechter theory and by adopting a binary 
scheme with accretion mass (Volonteri, Haardt & Madau 
2003). The evolution of gas and stars is then followed along 
the hierarchy with the following hypotheses: (a) stars can 
only form in haloes of mass > Msf{z); the evolution of 
Msf (z) will be discussed in the next Section; (b) in each halo 
the SF rate is proportional to the mass of cold gas whose 
gradual accretion is described by a numerically calibrated 
infall rate (Keres et al. 2005); (c) according to the so-called 
critical metallicity scenario (Schneider et al. 2002, 2006) low- 
mass stars with a Larson Initial Mass Function form when 
the gas metallicity Z > Zcr = 10~^^^ Zq; below Zcr massive 
Pop III stars form with a reference mass mpopiii = 2OOA'/0 . 
The chemical enrichment of the gas is followed by taking 
into account mass-dependent stellar evolutionary timescales 
(Lanfranchi & Matteucci 2007). We follow the chemical evo- 
lution of gas both in proto- Galactic haloes and in the GM 
by including a physically based description of the mechan- 
ical feedback due to supernova (SN) energy deposition (see 
SFS08 for details). Metals and gas are assumed to be instan- 
taneously and homogeneously mixed with the gas (implica- 
tions discussed in SFS08). The code is calibrated by simulta- 
neously reproducing the global properties of the MW (stel- 
lar/gas mass and metallicity) and the MDF of metal poor 
stars observed by Beers & Christlieb (2006) in the Galactic 
halo (for details see SSF07). 

Dwarf spheroidal candidates are selected among star 
forming haloes {Mh > M^f) which are likely to become 
satellites i.e. those corresponding to density fluctuations 
< 2a^ . Such a dynamical argument is based on N-body cos- 
mological simulations by Diemand, Madau & Moore (2005) . 
Therefore, at each redshift, we select dSph candidates from 
the merger tree among haloes with Adsf < Mh < M2cr; 
we then follow the isolated evolution (no further merging 
or accretion events) of virtual haloes with the same ini- 
tial conditions (dark matter/gas/stellar content, metallic- 
ity). Through this method we build a statistically significant 
dSph sample; however, it prevents us from making specific 
predictions on the actual number of dSph satellites. 



10 



* We adopt a ACDM cosmological model with h = 0.73, Qm = 
0.24, Ha = 0.72, n^h'^ = 0.02, n = 0.95 and erg = 0.74, consistent 
with WMAP5 data (Spergel et al. 2006); we assume the fit to the 
power spectrum proposed by Bardeen et al. (1986) and modified 
by Sugiyama (1995). 

t The quantity cr(M, z) represents the linear rms density fluctu- 
ation smoothed with a top-hat filter of mass M at redshift z. 
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Figure 1. Dark matter halo mass and circular velocity of se- 
lected dSph candidates as a function of their formation redshift z 
(points) for 10 realizations of the hierarchical merger tree. Differ- 
ent colours show the baryonic fraction at the formation epoch 
with respect to the cosmic value /c = 0.156: /f,//c > 0.5 (blue), 
0.1 < fblfc < 0.5 (green), fi,/fc < 0.1 (yellow). The lines show 
the evolution of Mgf(z) (solid), M^q^z) (dotted), the halo mass 
corresponding to 2cf peaks (dotted-long dashed), T^ir = 10'*K 
(short dashed line) and T^ir = 2 X lO'^K (long dashed line). 
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Figure 2. Total luminosity of selected dSph candidates as a func- 
tion of their average iron abundance (points) for 10 realizations 
of the hierarchical merger tree. As in Fig. 1 different colours show 
the baryonic fraction at the formation epoch with respect to the 
cosmic value. The points with error bars are observational data 
from Kirby et al. 2008. 



2.1 Feedback and star formation 

Radiative feedback processes are crucial to determine the 
minimum mass, Msf, and efficiency, /, , for star formation. 
The redshift evolution of the minimum mass of star-forming 
haloes, M^f, is determined by two distinct radiative feed- 
back processes. The first one has to do with the increase of 
the Jeans mass in progressively ionized cosmic regions; as 
a result, the infall of gas in haloes below a given circular 
velocity, w*, is quenched and such objects are gas-starved. 
The evolution of d* {z) depends on the details of the reioniza- 
tion history (Gnedin 2000, Schneider et al. 2008); we adopt 
the typical value v* = 30 kms~^ after the end of reionization 
Zrei = 6 I.e. we assume Msf{z) = M3o{z) when 2 < 6. Before 



© RAS, MNRAS 000, 1-5 



Ultra faint dwarfs: probing early cosmic star formation 3 




Figure 3. Left panel: comparison between tlie MDF of UFs observed by Kirby et al. (2007) (points) and tiie simulated one (histogram). 
Error bars are the Poissonian errors. The histogram is the averaged MDF over all UFs candidates (L < IO^Lq) in 10 realizations of the 
merger tree. The shaded area represents the la scatter among different realizations. Right panel comparison between the Sculptor MDF 
observed by Helmi et al. (2006) (points), and the simulated one (histogram). Error bars are the Poissonian errors. The histogram is the 
averaged MDF over all the Sculptor-like dSph candidates (IO^Lq < L < in 10 realizations of the merger tree. The shaded area 

represents the ztlcr Poissonian error. 



reionization a second type of feedback, related to the pho- 
todissociation of H2 molecules by the Lyman- Werner (LW) 
background photons, becomes important. As this species is 
the only available coolant for "minihaloes" {Tvir < 10"* K) 
SF is suppressed in these objects to an extent which depends 
on the intensity of the UV background (Haiman, Rees & 
Loeb 1996; Ciardi, Ferrara & Abel 2000; Kitayama et al. 
2000; Machacek, Bryan & Abel 2001). According to Dijk- 
stra et al. (2004) at 2 « 10 objects with Vc > 10 kms~^ 
{Tvir ~ 2000K) can self-shield and collapse; therefore we 
use this value as the minimum absolute threshold for star 
formation. During reionization {6 < z < 9) the interplay be- 
tween these two feedback types is quite complicated and no 
consensus is found on the evolution of Alsf (z) (see for exam- 
ple Fig. 25 of Ciardi & Ferrara 2005) . Instead of modeling in 
detail the build-up of LW and ionizing UV backgrounds, we 
interpolate between the low- and high-redshift behaviors and 
use an heuristic form for Msf {z) which leads to the suppres- 
sion of SF in gradually more massive objects (Fig. 1). Such 
parameterization of Msf{z) allows to correctly reproduce 
the observed iron- luminosity relation of dSphs (Fig. 2)"^. 
Fig. 1 also shows the halo masses of selected dSph candi- 
dates along with their initial baryonic fraction, /t,, with re- 
spect to the cosmic value fc = 0.b/^M = 0.156. We classify 
dSphs according to their initial baryonic content as gas-rich, 
h/ fc > 0.5, intermediate, 0.1 < ft/fc < 0.5, and gas-poor, 
fb/fc < 0.1, systems. A second important feature of the 
model is that we consider a mass- dependent SF efficiency. 
In minihalos, in fact, the ineffective cooling by H2 molecules 



t The total luminosity value is derived from the stellar mass con- 
tent as L = M* X (M/L), by assuming (M/L), = 1. 



limits the amount of gas than can be transformed into stars. 
Several authors (Madau, Ferrara & Rees 2001; Ricotti & 
Gnedin 2005; Okamoto, Gao & Theuns 2008) agree that in 
these systems the SF efficiency decreases cx T^i,.- A suitable 
form is then 



'2 X lO^K' 



(1) 



where e* is the value determined by matching the global 
properties of the MW. 



3 RESULTS 

The first test of the model is a comparison of the observed 
met allicity- luminosity relation with our predictions; the re- 
sults are reported in Fig. 2. The overall agreement is satisfac- 
tory, and the increasing trend of metallicity with luminosity 
well reproduced. The faint end of the relation, L < IVP Lq, 
is predominantly populated by minihaloes; above that lu- 
minosity H-cooling haloes dominate. The low-L segrega- 
tion of minihaloes naturally arises from their low masses 
and it is further extended towards very faint luminosities, 
L < 10''l/Q, by the reduced SF efficiency and, most impor- 
tantly, by their low initial baryonic fraction fb/fc < 0.1. 
Being UFs defined as systems with L < 10^ Lq, we conclude 
that all UFs are minihaloes but there are minihaloes that 
are not UFs. 

We have just pointed out that the initial baryonic frac- 
tion is a key factor in dSph galaxy evolution. What sets 
such initial value? Gas-rich systems are distributed through 



the entire mass (M^ 



10' 



Mq) and formation redshift 



{z = 2 — 10) ranges of selected dSph candidates (Fig. 1). 
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These dSphs are the present-day counterpart of predomi- 
nantly newly virialized objects accreting gas from the GM. 
Intermediate systems (0.1 < fb/fc < 0.5) originate from 
mixed merging of star- forming and starless progenitors; they 
are typically more massive than gas-rich haloes and form at 
lower redshifts. Their smaller baryonic content is the re- 
sult of shock-heating of the infalling gas during major merg- 
ing events (Cox et al. 2004) which stops accretion. Finally, 
gas-poor systems formed by merging of recently virialized, 
Mh > Msf progenitors. Since most of the diffuse gas is still 
accreting, shock heating quenches the infall early on. Note 
that a low baryonic content might also result from shock- 
stripping of the gas due to winds outflowing from nearby 
galaxies (Scannapieco, Ferrara & Broadhurst 2000). 

The faintest UFs are found to be gas-poor systems; this 
fact has two important implications. First, being the SN 
rate depressed in these systems by the scarce availability 
of gas, mechanical feedback has negligible effects, and they 
evolve as a closed-box. As essentially all metals are retained 
they have a relatively high Fe-abundance which is seen as 
an almost inverted Fe-L relation below L = lO^'^L©. Second 
UFs have extremely large M/L > 100 ratios, and the faintest 
ones among them reach such extreme values as M/L ~ 10'*. 

To get more insight in the nature of UFs it is instruc- 
tive to analyze their stellar MDF. In Fig. 3 (left panel) we 
show the MDF averaged over all UF candidates in 10 re- 
alizations of the merger tree, along with the available data 
(Kirby et al. 2008) . The shape of the distribution is relatively 
broad extending from [Fe/H]« —4 to ~ —1, in good agree- 
ment with observations. This is noticeably different from the 
Sculptor MDF shown in the right panel of the same Figure. 
SN feedback in Sculptor is much stronger than in UFs, as 
explained above; SF activity is then terminated shortly after 
100 Myr (see Fig. 3 of SFS08) whereas UFs can continue to 
quietly form stars up to 1 Gyr from their birth. These phys- 
ical differences are reflected in the wider UF MDF. In our 
model Sculptor-like dSphs-* are associated with gas-rich, H- 
cooling haloes, with Mh ~ 10* Mq, virializing at 2 ~ 7.5, in 
agreement with the findings in SFS08. However, in contrast 
with that study we now predict a small tail of [Fe/H]< —3 
stars. These are relics of the rare SF episodes occurred in 
some progenitor minihaloes (which were not considered in 
SFS08) at z > 7.5. Therefore, such a small number of very 
iron-poor stars in classical dSphs are expected to be char- 
acterized by the same abundance pattern of [Fe/H]< —3 
stars in UFs. The rest of the MDF, on the contrary, is built 
after the stellar population bulk assembling at lower red- 
shifts. A striking difference between the two distributions in 
Fig. 3 is the larger fraction (« 25%) of very metal poor stars 
([Fe/H]< -3) present in the UF MDF. The interpretation of 
this result is clear: as UFs form earlier the metallicity of the 
GM gas out of which they virialize is correspondingly lower 
but still high enough that low-mass stars can be produced 
according to the critical metallicity criterion. In principle, 
then, UFs are potentially powerful benchmarks to validate 
this criterion, which would exclude the presence of stars of 
total metallicity below Zcr, here fixed at Zcr ~ W~^'^Zq 
(see SFS08 for details). 



L W^ '^Lq, ([Fe/H]) fa -1.8, see also Fig. 2 
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Figure 4. Fraction = / fcM/^ of the potentially available 
cosmic baryon content turned into stars as a function of ■ The 
points are the average over all dSph candidates in 10 realizations 
of the merger tree. Error bars show the ilcr dispersion among 
different dSphs. The solid line is the handy fit. 



4 DISCUSSION 

UFs are the oldest and most dark matter [M/L > 100) 
dominated dSphs in the MW system with a total mass 
M = IQ^'^Mq and L = 10^"^Lq; they are found to be 
left-overs of H2-cooling minihaloes formed at z > 8.5, i.e. 
before reionization. Their MDF is broader (because of their 
more prolonged SF) and shifted towards lower [Fe/H] (as a 
result of a lower GM metallicity at the time of formation) 
with respect to classical dSphs. 

The SF in the faintest UFs is depressed by two factors: 
the limited availability of cold, star-forming gas due to inef- 
fective H2 cooling, and their specific formation mechanism 
resulting in gas-starved systems. The mass of stars formed 
at z = O normalized to the cosmic baryon content associated 
with a halo of mass Afh, — Mt/ f^Mh is shown in Fig. 4 
as a function of Mh ■ A handy fit to the curve is given by 

= exp{a + hx + cx^) (2) 

where a = -66, h = 5.97, c = -0.14 and x = ln{Mh/MQ). 
By inspecting that Figure the interplay between different 
factors affecting SF can be readily understood. Haloes above 
10® Mq convert about 10% of their potentially available 
fcMh baryonic mass into stars; this value is solely deter- 
mined by mechanical feedback. In haloes with 10* — 10® M© 
gas infall starts to be quenched by shock-heating during 
formation by mergers, with SN feedback playing a sub- 
dominant role. Therefore drops by about a factor 3 in 
this range. Finally, minihaloes Mh < 10* are affected by 
an additional suppression factor related to their decreasing 
SF efficiency (eq. 1) due to radiative feedback acting on the 
H2 chemical network. The larger scatter towards low Mh 
is hence due to the increasing number of physical processes 
influencing star formation. 

The recent work of Madau et al. (2008) attempts 
to determine for minihaloes by matching the luminos- 
ity function of MW satellites in the SDSS under the as- 
sumption that the formation of these systems stopped at 
Zrei = 11 due to rciouization photoheating. They find that 
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JF = (0, 0.0025, 0.02) for haloes with a total mass Mh = (< 
3.5 X W''Mq, (3.5 - 7) X IO^Mq, > 7 x IO'^Mq). In the same 
range of masses we find slightly higher values (Fig. 4); 
this difference is likely due to the assumption made here 
that Zrei = 6. It is however unclear what is the effect of an 
early, z > 11, formation on the median and shape of the 
minihaloes/UFs MDF; on the other hand our model cannot 
predict if the higher produces an excess of visible MW 
satellites. Clearly a more comprehensive study is required 
to answer these questions. 

Our conclusions support the "primordial scenario" for 
the origin of dSphs proposed by Bovill & Ricotti (2008). In 
such a scenario feedback effects of the type discussed here 
represent fundamental evolutionary ingredients. The picture 
is strengthen by the successful simultaneous match of the Fe- 
L relation and the MDF of UFs. Our model reproduces at 
the same time also the observed MDF of both the MW and 
of a prototypical classical dSph as Sculptor. 

The skeptic might wonder about the possible role of 
tidal stripping. In our picture UFs are, among all dSphs, 
those associated with the highest— it fluctuations (~ 2a) 
forming at the earliest possible epochs (Fig. 1). According 
to N-body cosmological sinmlations (Dicmarid, Madau & 
Moore 2005) this implies that UFs are most probably found 
at small Galactocentric radii, as indeed deduced from obser- 
vations; such data however must be interpreted with care as 
it might be biased due to the magnitude limit of the SDSS. 
The proximity to the MW may cause gas/stellar loss by 
tidal stripping from these satellites (Mayer et al. 2007). How- 
ever, it seems unlikely that UF can be the stripped remnant 
of classical dSph as the scaling of the luminosity-velocity 
dispersion with luminosity would be too steep to explain 
the observed trcud (Pcnarrubia, Navarro & McConnachic, 
2008). The success of our model also lends support to the 
conclusion that tidal stripping plays at most a minor role. 

A final caveat concerns the adopted heuristic assump- 
tion for Adsf{z). Physically Msf{z) is tightly related to the 
rcionization history of the MW environment. Although we 
have not attempted to model in detail the radiative feed- 
back processes determining the evolution of such quantity, 
which we defer to further study, it is conceivable that the 
grow of the LW background intensity will suppress the cool- 
ing and SF ability of progressively more massive minihaloes. 
Guided by this general argument we have then chosen an 
heuristic form of Ms/(z) which suitably accounts for this 
physical process as required by the data. A more physical 
interpretation of its shape derived from a detailed modeling 
of the LW background intensity growth during rcionization 
is nevertheless necessary; it may well unfold the complicated 
physics behind radiative feedback. 
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